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OBJECf IVE — Nonalcoholic fatty liver (NAFL) is thought to contribute to insulin resistance 
and its metabolic complications. However, some individuals with NAFL remain insulin sensitive. 
Mechanisms involved in the susceptibility to develop insulin resistance in humans with NAFL are 
largely unknown. We investigated circulating markers and mechanisms of a metabolically benign 
and malignant NAFL by applying a metabolomic approach. 

RESEARCH DESIGN AND METHODS— A total of 265 metabolites were analyzed be- 
fore and after a 9-month lifestyle intervention in plasma from 20 insulin-sensitive and 20 insulin- 
resistant subjects with NAFL. The relevant plasma metabolites were then tested for relationships 
with insulin sensitivity in 17 subjects without NAFL and in plasma from 29 subjects with liver 
tissue samples. 

RESULTS — The best separation of the insulin-sensitive from the insulin-resistant NAFL group 
was achieved by a metabolite pattern including the branched-chain amino acids leucine and 
isoleucine, ornithine, the acylcarnitmes C3:0-, C16:0-, and C18:0-carnitine, and lysophospha- 
tidylcholme (lyso-PC) C16:0 (area under the ROC curve, 0.77 [P = 0.00023] at baseline and 0.80 
[P = 0.000019] at follow-up). Among the individual metabolites, predominantly higher levels of 
lyso-PC C16:0, both at baseline (P = 0.0039) and at follow-up (P = 0.001), were found in the 
insulin-sensitive compared with the insulin-resistant subjects. In the non-NAFL groups, no 
differences in lyso-PC C16:0 levels were found between the insulin-sensitive and insulin-resistant 
subjects, and these relationships were replicated in plasma from subjects with liver tissue samples. 

CONCLUSIONS — From a plasma metabolomic pattern, particularly lyso-PCs are able to 
separate metabolically benign from malignant NAFL in humans and may highlight important 
pathways in the pathogenesis of fatty liver-induced insulin resistance. 
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The prevalence of nonalcoholic fatty 
liver (NAFL) is increasing worldwide 
and is now affecting >30% of adults 
and a considerable number of children in 
developed countries. NAFL represents a 
strong and independent predictor of type 
2 diabetes, the metabolic syndrome, and 
cardiovascular disease, findings that are 
probably largely attributable to dysregu- 
lated hepatic metabolic signaling, result- 
ing in insulin resistance, which is often 
found in fatty liver (1-8). In fact, in hu- 
mans, among several body fat compart- 
ments that are considered to play a major 
role in the pathogenesis of insulin resis- 
tance as visceral fat and intramyocellular 
fat, liver fat is most strongly associated 
with insulin resistance (9,10). However, 
a considerable amount of subjects with 
NAFL and without more advanced stages 
of fatty liver-associated diseases remain 
insulin sensitive (11). 

The following question remains: Can 
broad systematic blood screening help to 
find mechanisms and/or markers for the 
dissociation of NAFL from insulin resis- 
tance in humans? Furthermore, are these 
parameters also relevant in subjects with- 
out NAFL or can they even help to 
identify NAFL-associated insulin resis- 
tance? To address these questions in the 
current study, we used precise phenotyp- 
ing methods in humans that allowed 
careful quantification of body fat com- 
partments, liver fat content, and glucose 
and lipid metabolism and combined this 
setting with a targeted plasma metabolomic 
approach. 
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RESEARCH DESIGN AND 
METHODS 

Subjects 

Caucasians from the southern part of 
Germany participated in the ongoing 
Tubingen Lifestyle Intervention Program 
(12,13). Individuals were included in the 
study when they fulfilled at least one of 
the following criteria: a family history of 
type 2 diabetes, a BMI >27 kg/m 2 , and a 
previous diagnosis of impaired glucose 
tolerance and/or of gestational diabetes. 
They were considered healthy according 
to a physical examination and routine lab- 
oratory tests. The participants had no his- 
tory of liver disease and did not consume 
more than two alcoholic drinks per day. 
Serum aminotransferase levels were <2.5 
times the upper limit of normal. In a first 
approach, out of 330 subjects who met 
the aforementioned requirements, i.e., 
had measurements of body fat distribu- 
tion and liver fat content, using magnetic 
resonance techniques at baseline (12), we 
randomly selected 40 subjects who had 
NAFL for metabolomics analysis. We 
then measured the identified metabolo- 
mic parameters in 17 subjects who under- 
went the same phenotyping strategies and 
who did not have NAFL. Finally, we mea- 
sured the relevant metabolites in a third 
group of 29 subjects from whom liver tis- 
sue samples were available. Informed 
written consent was obtained from all 
participants after the nature and possible 
consequences of the studies were ex- 
plained and the local medical ethics com- 
mittee had approved the protocol. 

Lifestyle intervention 

The 57 subjects who underwent precise 
phenotyping participated in the 9-month 
lifestyle intervention. After the baseline 
measurements, individuals underwent di- 
etary counseling and had up to 10 sessions 
with a dietitian. Counseling was aimed to 
reduce body weight, intake of calories, and, 
particularly, intake of calories from fat and 
to increase intake of fibers. Individuals 
were asked to perform at least 3 h of 
moderate sports per week. Aerobic endur- 
ance exercise (e.g., walking or swimming) 
with an only moderate increase in the heart 
rate was encouraged (13). 

Total body fat, body fat 
distribution, and lean body mass 

Measurements of total body and visceral fat 
and lean body mass were performed by an 
axial Tl -weighted fast spin echo technique 
with a 1.5-T whole-body imager (14). 



Noninvasive measurement of liver 
fat content and intramyocellular 
lipids 

Liver fat content was noninvasively mea- 
sured by localized proton magnetic reso- 
nance ( HMR) spectroscopy as previously 
described (14). NAFL was defined as liver 
fat content >5.56% (15). Intramyocellular 
lipid content of the tibialis anterior muscle 
was determined as previously described (14). 

Oral glucose tolerance test and 
clinical chemical analyses 

The 57 individuals who underwent precise 
phenotyping and the intervention also 
underwent a 75-g oral glucose tolerance 
test (OGTT) at baseline and after 9 months 
of lifestyle intervention. Whole -body in- 
sulin sensitivity was calculated from glu- 
cose and insulin values during the OGTT as 
proposed by Matsuda and DeFronzo (16). 
In the 29 subj ects from whom liver samples 
were available, the homeostatic model as- 
sessment of insulin resistance (HOMA-IR) 
was calculated from fasting blood samples 
(glucose [mmol/L] • insulin [mU/L]/22.5). 
Blood glucose was determined using a bed- 
side glucose analyzer (YSI, Yellow Springs, 
CO). Plasma insulin was determined on 
an ADVIA Centaur XP and all other routine 
parameters on an ADVIA 1800 clinical 
chemistry system (Siemens Healthcare 
Systems, Erlangen, Germany). Serum plas- 
minogen activator inhibitor 1 was mea- 
sured by ELISA (Bender MedSystems), 
and fetuin-A was measured by an im- 
munoturbidimetic method (BioVendor 
Laboratory Medicine, Modreci, Czech 
Republic) (17). The metabolic syndrome 
was estimated based on the recommen- 
dation of the Executive Summary of the 
Third Report of the National Cholesterol 
Education Program (NCEP) Expert Panel 
on Detection, Evaluation, and Treatment 
of High Blood Cholesterol in Adults 
(Adult Treatment Panel III) (18). 

Liver samples 

The 29 Caucasians (62.8 ± 2.2 years of 
age; sex, 22 male/7 female; BMI, 25.5 ± 
0.85 kg/m 2 ) who underwent liver surgery 
were fasted overnight prior to the collec- 
tion of blood samples and liver biopsies. 
Subjects tested negative for viral hepatitis 
and had no liver cirrhosis. Liver samples 
were taken from normal, nondiseased tis- 
sue during surgery, immediately frozen in 
liquid nitrogen, and stored at — 80°C. 

Metabolite profiling 

At baseline and after 9 months of lifestyle 
intervention, 265 blood metabolites were 



analyzed by the targeted IDQ metabolo- 
mics platform from Biocrates (Innsbruck, 
Austria). This platform combines flow 
injection (acylcarnitines and glycero- 
phospholipids), liquid chromatographic 
(amino acids and bile acids), gas chro- 
matographic (free fatty acids), and mass 
spectrometric approaches. For the repli- 
cation of the findings on palmitoyl lyso- 
phospatidylcholine (lyso-PC C16:0) in 
the lifestyle intervention study, lyso-PC 
CI 6:0 levels were measured with another 
method, targeted ultraperformance liquid 
chromatography mass spectrometry, in 
the fasting blood samples of the subjects 
from whom liver samples were available. 

Statistical analyses 

First, 40 subjects with NAFL were ran- 
domly selected from a larger group solely 
based on their liver fat content at baseline, 
when they had NAFL (liver fat content 
>5.56%). Subjects were then divided 
into two groups based on the median in- 
sulin sensitivity, which was measured 
from the OGTT as proposed by Matsuda 
and DeFronzo (16) (insulin-sensitive and 
insulin-resistant NAFL groups). This sim- 
ple two-step approach was specifically 
chosen to not introduce a bias that may 
arise by an artificial matching process, 
e.g., when matching for age, sex, and or 
adiposity, because today it is not known 
how these parameters influence the rela- 
tionship of liver fat content with insulin 
sensitivity. The same approach was cho- 
sen for the 17 subjects without NAFL 
(insulin-sensitive and insulin-resistant 
non-NAFL groups). Pearson correlations 
were used to test for relationships be- 
tween continuous variables. Paired Stu- 
dent t tests were performed to compare 
data obtained from the same experimental 
group. When different groups were com- 
pared, two-sample Student t tests were ap- 
plied. In addition, repeated-measures 
multivariate AN OVA (MAN OVA) was 
used, which allowed us to test in a more 
general aspect for differences in the param- 
eters between the insulin-sensitive and 
insulin-resistant groups, as well as for 
changes of the parameters within the sub- 
jects and to perform group X time inter- 
action tests. For the evaluation of the mass 
spectrometric data, a signal-to-noise ratio 
of 3 was set as the limit of detection. If not 
otherwise stated, a significance threshold 
of a = 0.01 was used as a reasonable com- 
promise between type I (false-positive) 
and type II (false-negative) error, as previ- 
ously suggested (19). Statistical analysis 
was performed using Matlab version 
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7.9.0.529 (R2009b) (The MathWorks 
Inc.). 

To assess the discriminative potential 
of metabolite subsets in a multivariate 
manner, we applied a wrapper-based 
feature subset selection approach resort- 
ing to six different heuristic optimization 
techniques implemented in the EvA2 
framework and six classification algo- 
rithms from the WEKA package (20). 
We applied the feature subset selection 
procedures individually and additionally 
integrated them into an ensemble to dis- 
cover particularly robust metabolite sub- 
sets (21). Comparable approaches have 
demonstrated their usefulness in mass 
spectrometry studies (22). 

RESULTS 

Characteristics of the subjects with 
NAFL 

At baseline, the 20 insulin-sensitive and 
20 insulin-resistant subjects in the NAFL 
groups had almost an identical mean liver 
fat content, whereas the mean insulin 
sensitivity differed largely based on the 
allocation of the subjects into the groups 
(Table 1). There was a dominance of males 
(75%) in the insulin-sensitive group and 
of females (65%) in the insulin-resistant 
group as well as a significant difference 



in age between the groups. The insulin- 
resistant women were older than the 
insulin-sensitive women (54 ± 4 vs. 
46 ± 2 years, P = 0.06), and it is possible 
that more women in the insulin-resistant 
group were postmenopausal. General ad- 
iposity, as estimated from the BMI and 
precisely measured by magnetic reso- 
nance imaging, was lower in the insulin- 
sensitive group. No differences in plasma 
levels were detected for aspartate trans- 
aminase, alanine transaminase, high- 
sensitivity C-reactive protein (hs-CRP), 
fasting glucose, or 2-h glucose (Table 1). 
Also, no differences were found between 
the groups in the serum lipid profiles 
(Supplementary Table 1). Statistically 
significant higher fetuin-A levels, which 
may indicate elevated hepatic inflamma- 
tion, were found in the insulin-resistant 
group. During the lifestyle intervention, 
liver fat content decreased significantly 
and to a similar extent in both NAFL groups, 
and the difference in insulin sensitivity 
remained unaltered between the groups 
at follow-up (Table 1). In the repeated- 
measures MANOVA, the statistical rela- 
tionships at baseline and at follow-up, 
which were observed between and within 
the groups, were largely similar to the re- 
sults of the Student t tests (data not 
shown). 



Predictive effect of a metabolite 
pattern to separate insulin-sensitive 
from insulin-resistant subjects with 
NAFL 

After the initial evaluation of 265 plasma 
metabolite measurements, 85 metabolites 
had to be excluded based on low signal 
intensities (Supplementary Tables 2 and 
3). First, we hypothesized that a pattern 
of metabolites may be able to separate 
insulin-sensitive from insulin-resistant 
subjects with NAFL. The detected pattern 
included seven metabolites contributing 
to the separation of insulin-sensitive 
from insulin-resistant individuals with 
NAFL, both at baseline and at follow-up 
(Fig. 1A and D). We found that the sum of 
the levels of the metabolites of this pat- 
tern, which consisted of the branched- 
chain amino acids (BCAAs) leucine and 
isoleucine, ornithine, the acylcarnitines 
C3:0-, C16:0-, and C18:0-carnitine, and 
lyso-PC C16:0, was higher in the insulin- 
sensitive compared with the insulin- 
resistant NAFL groups at baseline and at 
follow-up (Fig. IB and E). More impor- 
tantly, using a naive Bayes classifier, this 
metabolite pattern at baseline separated 
insulin-sensitive from insulin-resistant 
subjects with a relatively high discrimina- 
tory power at baseline (area under the 
ROC curve, 0.77; positive predictive 



Table 1 — Selected characteristics of the 40 subjects with NAFL at baseline and after 9 months of lifestyle intervention 



Insulin-sensitive subjects Insulin-resistant subjects 



Characteristics 


Baseline 
(n = 20) 


Follow-up 


Baseline 
(n = 20) 


Follow-up 


P value at 
baseline 


P value at 
follow-up 


P-MANOVA 
group X time 


Insulin sensitivity (AU) 


14.7 ± 1.10 


18.4 


± 1.72 


7.3 ± 0.52 


7.9 ± 0.75 


<0.0001 


<0.0001 


0.75 


Liver fat MRS (%) 


8.6 ± 0.38 


5.8 


± 0.85* 


8.5 ± 0.51 


5.6 ± 0.65* 


0.72 


0.71 


0.53 


Sex (male/female) 


15/5 






7/13 




0.02 






Age (years) 
BMI (kg/m 2 ) 


52 ± 2 






44 ± 2 




0.016 






30.1 ± 0.76 


28.9 


± 0.93* 


32.4 ± 0.81 


31.8 ± 0.81* 


0.05 


0.02 


0.23 


Waist circumference (cm) 












0.61 


0.25 




Metabolic syndrome (%)# 


20 




22 


30 


50 


0.46 


0.07 




Total body fat M RT (kg) 


24.4 ± 1.82 


20.5 


± 2.27 


32.0 ± 1.50 


29.1 ± 2.06* 


0.002 


0.09 




Visceral fat MRT (kg) 


4.5 ± 0.41 


3.8 


± 0.54* 


3.7 ± 0.35 


3.6 ± 0.50* 


0.14 


0.83 


0.28 


LBM MRT (kg) 


61.3 ± 2.0 


59.3 


± 2.1 


58.6 ± 2.5 


60.9 ± 3.1 


0.40 


0.66 


0.77 


IMCL ti t> m li s anterior (AU) 


4.2 ± 0.6 


3.1 


± 0.3 


3.9 ± 0.4 


4.1 ± 0.4 


0.99 


0.04 


0.19 


AST (units/L) 


26.6 ± 1.46 


26.6 


± 3.44 


27.1 ± 2.58 


22.0 ± 1.51* 


0.89 


0.18 


0.16 


ALT (units/L) 


27.3 ± 2.33 


26.8 


± 4.58 


29.4 ± 3.30 


26.5 ± 2.99 


0.78 


0.95 


0.62 




















Fasting glucose (mmol/L) 


5.3 ± 0.15 


5.1 


± 0.15 


5.3 ± 0.12 


5.2 ± 0.12 


0.96 


0.86 


0.97 


2-h glucose (mmol/L) 


7.5 ± 0.37 


6.9 


± 0.50* 


7.3 ± 0.43 


6.4 ± 0.32* 


0.58 


0.51 




PAI-1 (ng/mL) 


4.5 ± 0.51 


4.6 


± 0.43 


5.5 ± 0.53 


6.5 ± 0.79 


0.10 


0.07 


0.25 


Fetuin-A (cig/mL) 


237 ± 9 






293 ± 14 




0.002 







Values are means ± SEM. ALT, alanine transaminase; AST, aspartate transaminase; AU, arbitrary units; IMCL, intramyocellular lipid; LBM, lean body mass; MRS, 
magnetic resonance spectroscopy; MRT, magnetic resonance tomography; PAI-1 , plasminogen activator inhibitor 1 . At follow-up, magnetic resonance measurements 
were only available in 30 subjects. #x 2 test. *P < 0.05 for change within each group. 
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Figure 1 — A plasma metabolite pattern predicts insulin sensitivity in NAFL. Plasma pattern of seven metabolites consistent with lyso-PC C16:0, 
ornithine, leucine, isoleucine, and C3:0-, C16:0-, and C18:0-carnitine for insulin-sensitive (blue) and insulin-resistant (red) subjects with NAFL at 
baseline (A) and after 9 months of lifestyle intervention (D). Each axis represents one metabolite, and all axes have their common origin in the center 
of the circle. The circle has radius 1 and intersects all axes at 1 (at the arrowhead) and —1 (in the lower half). The blue and red stars denote the 
individuals represented by the levels of the corresponding metabolites, projected to the axes after mean centering and scaling to unit variance. An 
individual's metabolite level was combined by equally weighted linear combination. The connecting lines between the stars meet at the center (i.e., 
mean) of all individuals of each group to indicate the separation of the groups in this projection. The closer the center of one of the groups is to I, the 
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value, 0.75; negative predictive value, 
0.88; sensitivity, 0.90; specificity, 0.70) 
(Fig. 1C). No such level of significance 
could be observed for any single metabo- 
lite. Interestingly, although many param- 
eters changed during the intervention, the 
metabolite pattern measured at follow-up 
also discriminated the insulin-sensitive 
from the insulin-resistant subjects at fol- 
low-up (positive predictive value, 0.83; 
negative predictive value, 0.77; sensitiv- 
ity, 0.75; specificity, 0.85) (Fig. IE and F). 

The metabolic fingerprints of 
insulin-sensitive and -resistant 
subjects with NAFL 

Next, we investigated which of the in- 
dividual metabolites were different be- 
tween the insulin-sensitive and the 
insulin-resistant subjects with NAFL. 
The plasma levels of total lyso-PCs 
reached the significance level P < 0.0 1 
at baseline (Fig. 2A). Within the metabo- 
lite group of lyso-PCs, lyso-PC C16:0 lev- 
els were different, with 19.4% lower 
plasma levels in insulin-resistant subjects 
(mean, 70.42 vs. 86.34 |xmol/L; P = 
0.0039) (Fig. IB). Because there was a 
sex difference between the NAFL groups, 
we also adjusted for sex using general lin- 
ear regression models, and the differences 
between the groups for total lyso-PCs and 
lyso-PC C16:0 were not largely affected 
(P = 0.008 and 0.015, respectively). Inter- 
estingly, the relationships were also not 
affected by the lifestyle intervention 
(e.g., for lyso-PC C16:0: insulin resistant, 
71.79 |jLmol/L vs. insulin sensitive, 
88.64 jjtmol/L), and the same differences 
in the metabolite parameters were found 
between the insulin-sensitive and the 
insulin-resistant NAFL groups after the 
intervention (Fig. 2C and D). 

Then, we investigated which metab- 
olites correlated with insulin sensitivity in 
all 40 subjects and in the 20 insulin- 
sensitive and in the 20 insulin-resistant 
NAFL groups analyzed separately. At 
baseline, the most robust correlations 
(significant relationships with insulin 
sensitivity were found in all subjects and 
in both NAFL groups analyzed sepa- 
rately) were present among the lyso-PCs 
(Supplementary Table 4). In all 40 sub- 
jects, for example, lyso-PC C16:0 
correlated not only positively with insulin 



sensitivity (r = 0.38, P = 0.016) but also 
negatively with circulating hs-CRP (r = 
-0.54, P = 0.0004), a marker of subclin- 
ical inflammation. During the lifestyle in- 
tervention, insulin sensitivity increased 
(P = 0.04) and hs-CRP levels decreased 
(P = 0.006), but lyso-PC CI6:0 levels 
did not change (P = 0.37). Interestingly, 
at follow-up, the relationship of lyso-PC 
C16:0 levels with insulin sensitivity was 
similarly strong (r=0.39,P = 0.014) com- 
pared with baseline; however, the rela- 
tionship of lyso-PC CI 6:0 levels with 
hs-CRP levels was no longer significant 
(r= -0.23, P = 0.16). 

Relationships of lyso-PC C16:0 levels 
with insulin sensitivity in subjects 
without NAFL 

We then tested whether the blood lyso- 
PC C16:0 levels were also different be- 
tween the respective non-NAFL groups 
(Supplementary Table 5). No difference 
in the lyso-PC C16:0 levels was found (in- 
sulin sensitive, 85.5 ± 3.3 |xmol/L; insu- 
lin resistant, 84.1 ± 6.7 |xmol/L; P = 
0.86). Power analyses revealed that we 
had a power of 0.97 to find the same dif- 
ference in lyso-PC C16:0 that we ob- 
served in the NAFL groups. 

Relationships of lyso-PC C16:0 levels 
with insulin sensitivity and hepatic 
inflammation in subjects with liver 
samples 

Finally, we investigated whether similar 
relationships were apparent in blood 
from subjects who underwent liver tissue 
sampling. Subjects in the upper (n = 7; 
triglyceride content, 4.8 ± 0.6%) and 
lower (n = 7; triglyceride content, 0.7 ± 
0.07%) quartiles of liver triglyceride con- 
tent were divided by the medians in fast- 
ing HOMA-IR. In the upper quartile, liver 
triglyceride content was not different 
(Supplementary Fig. 1A); however, the 
HOMA-IR was lower (Supplementary Fig. 
IB) and lyso-PC C16:0 levels were higher 
(Supplementary Fig. 1C) and a trend for 
less inflammation was found (Supplemen- 
tary Fig. 1D-F) in the insulin-sensitive 
compared with the insulin-resistant 
group. No difference in lyso-PC C16:0 lev- 
els was seen in the respective groups of 
the lower quartile (74.5 ± 2.4 vs. 64.8 ± 
2.4 |xmol/L, P = 0.56). 



CONCLUSIONS— Considering that 
among body fat and ectopic lipid com- 
partments, liver fat content correlates 
most strongly with insulin sensitivity 
(9,10), the fact that some subjects can ac- 
cumulate large amounts of fat within the 
liver yet remain insulin sensitive is an un- 
expected finding for many researchers. 
However, there is increasing data sup- 
porting the existence of such a phenotype 
(11,12,23), and the investigation of 
mechanisms determining this phenotype 
may help to advance our knowledge, spe- 
cifically about signaling and detoxifica- 
tion pathways that are involved in the 
regulation of insulin sensitivity in NAFL. 

In the current study, using a plasma 
metabolomic approach, we first found 
that a metabolite pattern including the 
BCAAs leucine and isoleucine, ornithine, 
the acylcarnitines C3:0-, C16:0-, and C18:0- 
carnitine, and lyso-PC CI 6:0 was able to 
separate the insulin-sensitive from the 
insulin-resistant NAFL groups with a rela- 
tively high discriminative power. Interest- 
ingly, the same metabolomic pattern 
separated the NAFL groups after a lifestyle 
intervention, after liver fat content and 
insulin sensitivity had changed, yet similar 
relationships for liver fat and insulin sensi- 
tivity were still apparent between the 
groups. When we then studied the individ- 
ual metabolites comprising this metabolo- 
mic pattern, we found lyso-PCs to most 
strongly correlate with insulin sensitivity 
in subjects with NAFL. These findings 
support that a metabolomic pattern, partic- 
ularly one including lyso-PC C16:0, may 
be a circulating marker of insulin sensitivity 
in NAFL. 

The metabolomic pattern that we 
have identified to separate the NAFL 
groups includes BCAAs. They were 
shown by interesting studies from New- 
gard et al, Gerszten, and others (19,24- 
26) to be important for the regulation of 
glucose metabolism. In addition, Perse- 
ghin and colleagues (27) showed that in 
first-degree relatives of type 2 diabetic in- 
dividuals, not only glucose and fatty acid 
metabolism but also leucine metabolism 
showed resistance to insulin. In our pres- 
ent study, we also detected correlations of 
leucine and isoleucine with insulin sensi- 
tivity at baseline in our 40 subjects with 
NAFL. However, and in contrast to what 



higher the levels of the metabolites in that group compared with the comparator group. Differences in this plasma metabolite pattern between insulin- 
sensitive and insulin-resistant subjects with NAFL at baseline (B) and after 9 months of lifestyle intervention (E ); depicted P for differences between the 
groups in two-sample Student t tests. P < 0.0001 for a group, P = 0.26 for time, and P = 0.55 for a group X time effect in the repeated-measures 
multivariate ANOVA test. Discriminating power of the plasma metabolite pattern to differentiate insulin-sensitive from insulin-resistant subjects with 
NAFL at baseline (C) and after 9 months of lifestyle intervention (?) (area under the ROC curve). AU, arbitrary units; AUC, area under the curve. 
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Figure 2 — Plasma lysophosphatidylcholine (lyso-PC) levels in insulin-sensitive and insulin- 
resistant subjects with NAFL. Plasma levels of total lyso-PCs (A) and lyso-PC CI 6:0 (B) of the 
insulin-sensitive group with NAFL and the insulin-resistant group with NAFL before the 9 months 
of lifestyle intervention (baseline). Plasma levels of total lyso-PCs (C) and lyso-PC C16:0 (D) 
after the 9 months of lifestyle intervention of the insulin-sensitive group with NAFL and the 
insulin-resistant group with NAFL; depicted P for differences between the groups in two-sample 
Student t tests. P = 0.001 for a group, P = 0.27 for a time, and P = 0.44 for a group X time effect 
(total lyso-PCs) (A and C) and P = 0.0009 for a group, P = 0.35 for a time, and P = 0. 65 for a group 
X time effect (palmitoyl lyso-PCs) (B and D) in the repeated-measures multivariate ANOVA test. 



we had expected, these relationships were 
positive. Interestingly, in the smaller 
group of 17 subjects without NAFL, these 
relationships were negative (both r < 
-0.42 and both P < 0.09). In addition, 
we detected a strong interaction effect of 
these circulating BCAAs with NAFL (both 
P < 0.006) but not with BMI, to deter- 
mine insulin sensitivity. Although we 
have no explanation for this novel ob- 
servation, and cannot rule out that the 
relatively small sample size may have 
brought about this finding, it may be of 



importance for further studies on the role 
of BCAAs in the regulation of insulin sen- 
sitivity, particularly in the aspect of 
NAFL. 

What is known about the role of lyso- 
PCs in metabolism? Lyso-PCs in plasma 
are primarily generated from PCs of lip- 
oproteins or from membrane-derived PCs 
by the action of secretory or lipoprotein- 
bound phospholipase A 2 . Another 
source of plasma lyso-PC levels is the pro- 
duction from PCs by endothelial lipase 
as well as from HDL or oxidized LDL 



by lecithin-cholesterol acyltransferase, 
which is secreted from the liver. On the 
other hand, lyso-PCs can also be reacyla- 
ted to PC by the action of lyso-PC acyl- 
transferase. Interestingly, most recently, a 
study in mice identified palmitoyl-, 
stearoyl-, and oleoyl-lyso-PCs to be lower 
in the serum of NASH animals compared 
with animals with simple steatosis and 
suggested that high lyso-PC acyltransfer- 
ase activity is involved in this process 
(28). By the action of lysophospholipase 
D/autotoxin plasma, lyso-PC can also be 
transformed to lysophosphatidic acid 
and sphingosine-1 -phosphate, important 
extracellular signaling lipids with multi- 
ple biological functions (29). Finally, the 
bioavailability of choline, which is re- 
quired for hepatic PC genesis, plays an 
important role in the pathophysiology 
of NAFL and its progression, and 
is regulated by the gut microbiome (30), 
may determine the plasma lyso-PC levels. 
Interestingly, recently lyso-PC C16:0 was 
found to enhance glucose uptake in an 
insulin-independent and protein kinase 
C-8-dependent manner in adipocytes 
(31). Lhis finding was confirmed in vivo 
by glucose-lowering effects of lyso-PC 
C16:0 in type 1 and type 2 diabetes 
mouse models (31). Whether these ef- 
fects and putatively beneficial effects of 
lyso-PC C16:0 on glucose uptake in other 
tissues translate into a positive effect on 
whole-body glucose uptake needs to be 
investigated in future studies. 

In contrast to glucose metabolism, 
the role of lyso-PCs in the regulation of 
inflammation has been studied more ex- 
tensively. It has been shown that lyso-PC 
CI 6:0 is able to inhibit reactive oxygen 
species production in stimulated poly- 
morphonuclear leukocytes (32,33). As 
putative mechanisms of action, lyso-PCs 
activate G protein-coupled receptors 
(GPR4, G2A, and GPR119) that are in- 
volved, for example, in migration and apo- 
ptosis of immune cells, osteoclastogenesis 
and angiogenesis, and insulin secretion 
from (3-cells (29). Importantly, lyso-PCs 
were recently found to enhance the sup- 
pressive function of human CD4 + CD25 + 
regulatory L cells (L regs ) (34). L regs secrete 
anti-inflammatory signals, inhibit macro- 
phage migration, and induce noninflamma- 
tory M2-like macrophage differentiation 
(35). In addition, defects in L reg function 
are considered to play a role not only in 
autoimmune diseases, such as autoimmune 
hepatitis (36), but also in other causes of 
hepatic inflammation (37,38). In agreement 
with an anti-inflammatory property of 
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lyso-PCs, we detected strong negative rela- 
tionships of total lyso-PCs and lyso-PC 
C16:0 with hs-CRP levels in our subjects 
with NAFL. 

We then asked whether lyso-PC C16:0 
levels also differ between insulin-sensitive 
and insulin-resistant subjects without 
NAFL. No such relationship was found. 
Importantly, the insulin-sensitive and 
insulin-resistant subjects in the respective 
NAFL and non-NAFL groups had very 
similar insulin sensitivity (Fig. 3). This 
indicates that lyso-PC C16:0 levels may 
be a good marker of subclinical inflam- 
mation and insulin resistance and may 
putatively protect from inflammatory pro- 
cesses, specifically in NAFL. In support of 
this hypothesis, we found lyso-PC C16:0 to 
be negatively associated with insulin resis- 
tance in subjects with high, but not low, 
liver tissue lipid content. Certainly, because 
of the very small sample size of this group, 
we cannot draw definite conclusions from 
this observation. Importantly, the replica- 
tion was not only done in a separate 
population than the initial screening pop- 
ulation but another method of measure- 
ment of lyso-PC C16:0 levels was also 
applied. Interestingly, in the large Relation- 
ship Between Insulin Sensitivity and Car- 
diovascular Risk (RISC) study, lyso-PCs 
were also found to negatively correlate 
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with estimated hepatic insulin resistance 
(39). Also in agreement with our present 
data, the correlations were stronger in the 
subjects with NAFL compared with the 
subjects without hepatic steatosis (A.G., 
personal communication). Moreover, and 
in agreement with our hypothesis, our liver 
histology data indicate that hepatic inflam- 
mation is present when plasma lyso-PC 
C16:0 levels are low. 

In summary, from a panel of 180 
metabolites, the lyso-PC concentrations 
differed most strongly between insulin- 
sensitive and insulin-resistant subjects 
with NAFL. These relationships, as well as 
the correlations with insulin sensitivity, 
were still apparent after 9 months of 
lifestyle intervention. Interestingly, lyso- 
PC levels did not change during this 
intervention, indicating that they are ge- 
netically determined. The fact that the 
correlations of lyso-PC levels with insulin 
sensitivity and with systemic and hepatic 
inflammation were less strong or absent 
when liver fat content was low, and the 
knowledge about the signaling properties 
of lyso-PCs, supports the hypothesis that a 
genetically determined high lyso-PC avail- 
ability may protect specifically from fatty 
acid-induced insulin resistance and he- 
patic inflammation when NAFL is present. 
Alternatively, elevated circulating lyso-PCs 
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Figure 3 — Relationships of liver fat content with insulin sensitivity at the baseline visit. The insert 
depicts the mean (SE) values for the respective groups and the level of statistical significance. 
*P < 0.05 for differences between the insulin-resistant and insulin-sensitive groups. P = 0.20 for 
differences between the insulin-sensitive and P = 0. 1 8 for differences between the insulin-resistant 
groups. #P < 0.05 for differences between the insulin-resistant NAFL group and all other groups. 
AU, arbitrary units. 



may not be functionally active in this pro- 
cess but represent activated pathways in the 
synthesis or metabolism of phospholipids 
and phosphatidylcholine and, thus, in the 
generation of hepatic endoplasmic reticu- 
lum stress (40) (Supplementary Fig. 2). 

In conclusion, in the future, a metab- 
olomic fingerprint, most robust in the class 
of lyso-PCs, may be able to differentiate 
insulin-sensitive from insulin-resistant sub- 
jects with NAFL. Furthermore, as low lyso- 
PC levels appear to specifically be a marker 
of NAFL-associated insulin resistance, 
these findings may highlight novel and 
interesting pathways for the studies on the 
pathogenesis of insulin resistance in NAFL. 
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